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A method is proposed for measuring the electron density at known points in the outer 
ionosphere, by the use of vlf receiving equipment in an artificial satellite, in conjunction 
with a vlf transmitter on the ground. The transmitter would radiate continuous waves, 
which would be propagated through the ionosphere in the "whistler" mode. The basis 
of the method is a measurement of the local wave admittance of the medium, by comparison 
of the signals received on an electric dipole and on a loop. 

A further proposal is made for an integrated vlf satellite experiment, in which several 
different types of observation would be made simultaneously. 



1 . Introduction 

Until recently, most of our knowledge of (he 
ionosphere bad boon obtained by indirect methods, 
such as vertical-incidence soundings with radio waxes. 
Now, however, the ionosphere can he investigated 
directly, with equipment borne aloft in rockets or 
satellites. These (wo vehicles lend themselves to 
making complementary studies: a rocket, which goes 
straight up and then falls back (o earth, is better for 
studies of the variations of the properties of the 
ionosphere 4 with height at a particular time and 
place, while a satellite, which may stay in orbit for 
a long period, is better for synoptic studies of slow 
variations with time and latitude at particular 
heights. Rocket and satellite experiments are 
especially valuable as a source of information about 
the outer ionosphere, above the peak of the F2 layer, 
which is beyond the reach of conventional vertical 
soundings. 

In studying the ionosphere, the number density 
of free electrons is a quantity of major interest. 
Already there are several ways in which the influence 
of the ionospheric electrons has been observed on 
transmissions from the artificial satellites; examples 
are the Faraday rotation of the plane of polarization, 
and the difference in refraction at disparate frequen- 
cies. Such observations, however, all measure 
integrated effects of the distribution of electrons 
along the entire path from the satellite to the 
observer, rather than their density at any one point. 
But the special advantage of a satellite is that it 
provides the opportunity for making local measure- 
ments at known points in the ionosphere, and so the 
question arises of how to equip it to measure local 
electron density. This paper sets out a proposal for 
one new method. Elsewhere, other methods have 
been suggested (Hoffmann [l]). 1 Since few have been 
put to practical test, it is too early to assess their 
relative merits. 



i Figures in brackets indicate the literature references at the end of this paper. 



The basis of the proposed method is the measure- 
ment of the local wave admittance of the ionospheric 
medium. To understand this concept, consider first 
a linearly-polarized plane wave in free space. The 
electric (E) and magnetic (H) vectors, which com- 
prise the field of the wave, lie at right-angles to each 
other in the plane of the wavefront (fig. 1). At any 
fixed point, the variations of the two vectors are 
always in phase, so that their ratio E/HisB, constant. 
This constant, which has the dimensions of a re- 
sistance, is called the wave impedance of free space 
(AV); its value is about 'M7 ohms. Its physical 
significance may he appreciated from the fact that if 
the wave were to impinge directly upon a uniform 
sheet of lossy material, that had this resistance per 
unit square area of surface, then it would he ab- 
sorbed without reflection. Thus such a sheet would 
be "matched" to free space. 2 
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Figure 1. A linearly-polarized plane wave in free space. 



2 Strictly, for the absorption to be complete the resistive sheet would have to 
be backed by a perfectly reflecting sheet, at a distance of a quarter-wavelength. 
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Within a uniform dielectric medium, the ratio of 
E to H is different from its value in free space. It is 
given by 

E/H=B=RJn (1) 

where n is the refractive index. In words, the wave 
impedance of the medium is inversely proportional 
to its refractive index. Henceforth, therefore, it will 
be more convenient to talk in terms of the wave 
admittance, which is the reciprocal of the impedance, 
and so is directly proportional to the refractive 
index : 

H/E=A=vA . (2) 

Evidently, if the wave admittance of the medium 
could be measured in some way, this measurement 
would yield the value of the local refractive index. 

In applying this result to the measurement of local 
electron density in the ionosphere, the following 
experimental arrangement is envisaged: There is to 
be a transmitter on the ground, radiating continuous 
waves of constant amplitude. The waves travel 
upward into the ionosphere, where an orbiting 
satellite is equipped to receive them. The satellite 
receives the waves on two separate antennas, a loon 
for the magnetic field and a dipole for the electric 
field. Information about the strength of the two 
fields is telemetered back to ground. At the ground 
station, the ratio of the two field strengths is com- 
puted, giving the wave admittance of the ionosphere 
in the immediate neighborhood of the satellite. By 
comparing the observed wave admittance with the 
value to be expected for free space, the local refrac- 
tive index is obtained. Now the refractive index is 
related to the 4 electron density by the expressions of 
the magneto-ionic theory. Hence, given the local 
refractive index, it is possible to work backwards 
and arrive finally at the value of the local electron 
density. These are the essential features of the 
proposed experiment, though there are many prob- 
lems of detail. 

The first problem to consider is the choice of the 
frequency of the transmitted waves. In principle, 
the experiment would work at any frequency. But 
if measurements are to be made above the maximum 
of the F2 layer, then the choice is restricted to those 
frequencies that are capable of passing through this 
layer: either the frequency must exceed the "critical 
frequency" of the layer, which varies roughly in the 
range 5 to 15 Mc, or it must be less than the gyro- 
frequency (about 1.5 Mc), in which case the waves 
can penetrate the layer in the "whistler" (ordinary) 
mode. It is preferable to choose low frequencies 
and the "whistler" mode, on the following two 
counts : 

First, at low frequencies only this mode is propa- 
gated, whereas at high frequencies both the ordinary 
and extraordinary mode would be propagated; 
their mutual interference would produce a compli- 
cated pattern of field in the ionosphere. 

Second, at low frequencies the waves are affected 
very strongly by the ionospheric electrons, so that 
the refractive index differs greatly from unity; high 



frequencies are affected much less. To illustrate 
this point, suppose that at some level above the F2 
layer the electron density is such as corresponds to 
a plasma frequency of 1.5 Mc (see sec. 2.1.). For 
the transmitted frequencies, take 15 kc as a typical 
low value and 15 Mc as a typical high value. The 
refractive indices for these two frequencies are then 
as given in table 1 below: 

Table 1. Refractive indices of the ionosphere, for a plasma 
frequency of 1.5 Mc 



Wave frequency 


Refractive index 


15 kc 


12. 5 (whistler mode). 


15 Mc 


0. 99 (ignoring magnetic field). 







In this example, the change of the refractive index 
from unity is more than a thousand times greater at 
the lower frequency. There can be no doubt, there- 
fore, that low frequencies are preferable. However, 
since the "whistler" mode requires the presence of 
a magnetic field to support its propagation, the use 
of this mode does confine the experiment to the 
domain of the earth's magnetic field. 

Having chosen the frequency and mode of propa- 
gation of the transmitted waves, the remaining 
problems arise from the fact that, at low frequencies, 
the ionosphere does not behave at all like a simple 
dielectric, so that the whistler mode has quite a 
complicated field structure. The outstanding ques- 
tions are: (a) What is the structure of the wave 
field in the whistler mode, and how is the wave ad- 
mittance defined for this mode? (b) How should the 
satellite be equipped to explore the structure of 
the field, and measure the wave admittance? (c) 
What sort of records would be obtained at the 
ground? (d) How should the records be interpreted? 
These questions are discussed successively in sec- 
tions 2 to 5 of this paper. Section 6 describes a 
proposal for an integrated vlf satellite experiment, 
in which several different types of observation 
would be made simultaneously, so as to obtain the 
fullest information possible from the vlf transmis- 
sions; it includes an independent method for ob- 
taining another measure of the local electron density. 
Finally, section 7 reviews the proposals and ex- 
amines some possible sources of error. 

2. Structure of Waves in the Whistler Mode 

2.1. Introduction 

In the proposed experiment, a transmitter on the 
ground radiates continuous waves at a very low 
fixed frequency. Some of these waves reach the j 
lower boundary of the ionosphere, where their 
energy is divided between the two magneto-ionic 
modes. The ordinary mode is reflected, with some 
absorption. The extraordinary mode also is 
partly reflected and partly absorbed, but in addition 
it experiences a mode coupling effect around the 
ordinary reflection level. The coupling is brought 
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aboul by collisions, and by the rapid variation of 
the electron density with height in the ionosphere at 
(his low level. Its effect is to transfer some energy 
from the extraordinary mode through to that 
hianch of the ordinary that is propagated freely at 
large electron densities; this branch is the mode of 
propagation of "whistlers." Judging from measure- 
ments on the "whistlers" (Storey [2]), energy is 
coupled into the ionosphere over an area around 
the transmitter that is about 1,000 km in radius. 
Thus, there is a vlf field in the ionosphere above 
this whole area, due to waves that have penetrated 
the lower boundary and traveled upwards in the 
"whistler" mode, together perhaps with some down- 
coming waves, that would arise from partial reflec- 
tions at any steep gradients of electron density. 
In the upper atmosphere, where the medium is 
probably fairly uniform, there is likely to be just a 
single upgoing plane wave. Moreover, owing to 
the strong refraction at the base of the ionosphere, 
the wave is likely to be traveling more or less radially 
outwards from the earth (Storey [2]). The structure 
of the field of this wave must be considered now. 

In describing the field structure, the rectangular 
coordinate system that is shown in figure 2 will be 
adopted. The coordinate axes are labeled 1, 2, and 
3, and in (hat order form a right-handed set: axis 3 
lies along the direction of propagation, that is to say, 
the wave normal; axis 1 is directed so thai the geo- 
magnetic field vector II lies in the 1 :> plane, and 
has a positive 1 -component : axes 1 and 2 define the 
plane of the wavefront. The various components of 
the vectors H and E in the field of the wave will be 
indicated thus: H=(fli, H 2 , H z ) and E (E u E 2j E z ). 
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Figure 2. The coordinate system. 



The term "structure" is used here to mean the 
relative values of the various axial components of the 
wave field at a point, and how these components vary 
from point to point and from one instant of time to 
another. The structure, in this sense, is independent 
of the amplitude of the wave, because the medium is 



linear. For determining the electron density, the 
absolute amplitudes of the field components need not 
be known, and the only requirement is that they be 
large enough to overwhelm external atmospheric 
uoise and circuit noise in the receivers. Neverthe- 
less, the absolute amplitudes would be measured in 
the course of t he experiment, and they are of interest 
for other reasons. 

The structure of a, plane wave of constant fre- 
quency can be described fully by the following 
quantities : 

(a) Refractive index. If there is no absorption, 
all components of the wave vary in time and space as 

exp {jo)(t — fiX 3 /c)} 

where n is the refractive index, w(^27r/) is the angular 
frequency of the wave, t is the time, c is the speed of 
light in free space, and x z is the component of dis- 
tance along the wave normal. 

(b) Polarization. The polarization is specified by 
the ratios of the components of II among themselves 
(//, :// :< : // :v ), and of the components of E among 
themselves (h\: IC>: E z ). 

(c) Wavi admittance. The wave admittance is the 
ratio of one component of IT to one of E. Here use 
will be made of the admittance looking in the direc- 
tion of propagation, which is defined as //, K { . 

These quantities are functions of the frequency of 
the wave, its direction of propagation relative to 
the earth's magnetic field, and the local properties 
of the medium. Expressions for their values are 
provided by the magneto-ionic theory; the expres- 
sions quoted here are taken, with some change of 
notation, from the work of Booker [3]. Rationalized 
units are used, with the following definitions: 

( bnstants: 

e, charge of electron (a negative quantity) 
///, mass of electron 
e n , electric permittivity of free space 
juo, magnetic permittivity of free space. 

\ ariables: 

/, frequency of wave 

6, angle between wave normal and earth's magnetic 

field 
A r , number density of electrons 
H , strength of earth's magnetic field. 

Characteristic frequencies: 

(„ is the electron plasma frequency. 

f H is the electron gyrofrequency. 
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Dimensionless xjarameters: 

Y=a) H /u=f H /f 

Y L =Y cos0 
Y T = Y sin0. 

With the aid of this notation, the refractive index, 
polarization, and wave admittance are discussed in 
sections 2.2, 2.3, and 2.4, respectively. The influence 
of collisions is ignored, because they are infrequent 
in the outer ionosphere. 

To illustrate the results of the theory, and to 
examine the adequacy of any approximations, nu- 
merical calculations will be made for the following 
particular set of conditions : 



/=15.5kc 
(9=160° 



/ p =1.55Mc 
/^=1.06Mc. 



The parameters of the magneto-ionic theory then 
have the values 



F=68.7 



Fi=-64.5 
F r =23.5. 



Such conditions might apply at a height of about 1000 
km above the vlf radio transmitter NSS (geomag- 
netic latitude 50° N). The wave is assumed to be 
traveling straight out from the earth in the plane of 
the magnetic meridian, so that axis 3 of the coordi- 
nate system is vertical, while the other two axes are 
horizontal; axis 1 points towards magnetic north, and 
axis 2 towards magnetic east. Note that in this 
instance the angle 6 is greater than 90°, since the 
earth's magnetic field is directed downwards in the 
northern hemisphere. 

2.2. Refractive Index 

When collisions are ignored, the general expression 
for the refractive index becomes 



ic 2 =l- 



2X(\-X) 



2(\-X)-Yl± Vr£+4F£(1-X) 2 



(3) 



in which the positive sign corresponds to the ordinary 
mode, and the negative to the extraordinary. At 
wave frequencies that are low compared to both the 
local plasma frequency and gyrofrequency, the re- 
fractive index for the "whistler" (ordinary) mode is 
usually given quite accurately by the quasi-longi- 
tudinal (QL) approximation: 



M 2 =l- 



X 



(4) 



The condition for the approximation to hold is that 

l(x-i)iy»m. (5) 



This equality is well satisfied in the example chose]], 
where the error of the QL approximation to the re- 
fractive index is less than 0.1 percent. 
The further approximation 



X 



SI 



\Y L \ ff B \ cos e\ 



(6) 



holds if X»\Y L \»1. In the example, this 
approximation gives jjl with an error of about 1 
percent. 

2.3. Polarization 

Consider first the polarization in the plane of the 
wavefront. The ratios of the transverse components 
of H and E among themselves are given by the 
"polarization ratio" 



R=-H 2 IH 1 =E l /E 2 



where 



and 



= -jk{l±^l+k~ 



k=: 



-X 



1 F| = l y sin^ 
= 2Y L ~2 cosfl 



(7) 
(8) 

(9) 
(10) 



and where the positive and negative signs cor- 
respond to the ordinary and extraordinary modes 
respectively. 

Now the "whistler" mode is the branch of the 
ordinary mode on which X is greater than unity, so 
that the denominator in (9) is negative; hence k is 
negative or positive according to whether d is less or 
greater than 90°. Also, at vlf the quasi-longitudinal 
condition (5) is satisfied, and this condition implies 
that \k\ <C 1. With these restrictions on the 
value of k, the polarization ratio for the whistler 
mode is approximately 



R^+j{l+\k\} 



(11) 



where the sign is positive or negative according as 6 
is less or greater than 90°, that is, according as the 
direction of propagation makes an acute or obtuse 
angle with the direction of the earth's magnetic field. 
In general, the polarization is elliptical, with tin 4 
axes of the ellipse parallel to axes 1 and 2 of the 
coordinate system; this conclusion follows from the 
fact that R is purely imaginary, so that H l and H 2 
are in phase quadrature. The vector sum traces 
out an ellipse with its minor axis parallel to axis 
1 . The sum of E x and E 2 traces out an ellipse of the 
same shape, but in this case the minor axis is parallel 
to axis 2. 
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The axial ratios of these ellipses are very close to 
unity; in the example, they differ from unity by less 
than 0.1 percent. So the polarization is almost 
circular, as it always is under QL conditions. 

The sense of rotation around the circle is specified 
by the sign of the imaginary value of R, a positive sign 
meaning right-handed rotation. For the whistler 
mode, the rotation is right-handed when 6 is less 
than 90°, and left-handed when 6 is more than 90°. 
Of course, the senses of rotation for the H and E 
vectors are the same. 

Now consider the longitudinal field components, 
the components along axis 3. The H vector has no 
such component; indeed, radio waves in a plasma 
have the general property that their magnetic field 
is always wholly transverse to the direction of prop- 
agation. Their electric field, however, in general 
does have a longitudinal component, which is 
given by 



E,/E 2 =j 



;(M 2 -1)F T 



1-X 



A!^!-i)U-f) 



(12) 



(13) 



The approximation holds under 
Furthermore, if X » 1 and \Y L \ » 



QL conditions. 



E,/E 2 ^- 



J \Yl\ 



= — il 



;|tan0|. 



To tl 



ie same approxima! ion 



Es/E^ — tsmd. 



(14) 



(15) 



These expressions mean that 2? 3 is in phase with 
E h and of such a magnitude and sign that the sum 
of these two components is roughly perpendicular 
to the geomagnetic field; in the numerical example, 
the angle between the sum and the perpendicular is 
about 3/4°. Now E 2 also lies in the plane per- 
pendicular to the field, so that the total E vector- 
traces out an ellipse in this plane, with the minor 
axis parallel to axis 2 of the coordinate system. The 
explanation of this result is that the electrons in the 
ionosphere are able to move much more rapidly in 
the direction of the geomagnetic field than at right 
angles to it, and at low frequencies they have 
sufficient time during a period of oscillation to move 
and neutralize any component of the electric field 
in this direction. Thus the result is not true in 
general, but only at low frequencies. 



2.4. Wave Admittance 

The admittance looking in the direction of propa- 
gation is defined as 



A=H 2 /E, .HJE 2 . 



(16) 



1 1 is given by 



A=nA 



(17) 



where /* is the refractive index, given by (3), and A 
is the admittance of :" 



tree space 



A = (cq/mo) *= (377 ohms) " 



(18) 



These expressions show that the held components 
H 2 and E x are in phase. From this and previous 
results it follows that, in the plane of the wavefront, 
the rotating magnetic vector leads the projection of 
the electric vector by 90° when there is propagation 
along the earth's magnetic field (0<9O°), and lags 
it bv 90° when the propagation is against the field 
(0>9O°). 

2.5. Summary of Approximations 

To summarize, the structure of w r aves in the 
whistler mode can be represented quite well by the 
following approximate expressions: 



For the refractive index, 



X 



ii" =i 



IFJ-1 



For t he polarization, 



E, = ±jE 2 
Ki= —E x tanfl 



(19) 



(20) 
(21) 
(22) 



where the sign in (20) and (21) is positive or negative 
according as 6 is less or greater than 90°. The 
assumptions now are thai the H vector traces out 
an exact circle in the plane of the wavefront, and 
the E vector an ellipse in the plane exactly per- 
pendicular to the earth's magnetic field, such that 
the projection of the ellipse onto the plane of the 
wavefront is a circle also. This state of polarization 
is illustrated in figure 3; here 9 has been taken as 
acute, for the sake of clarity, whereas this angle is 
obtuse in the worked example. 
Finally, for the wave admittance, 



A=H 3 /E l =pA 



(23) 



with n given by (19). 

These approximations are used throughout the 
rest of this paper. 

3. Details of the Satellite 
3.1. General Description 

The satellite would be equipped with two separate 
antennas for receiving the vlf waves. One won Id 
be a loop, responding to the magnetic field of the 
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Figure 3. The polarization of waves in the whistler mode. 



wave, while the other would be a dipole, for the 
electric field. The two antennas would be mounted 
coaxialry, that is to say, the plane of the loop would 
be perpendicular to the axis of the dipole. With 
this arrangement, the satellite would receive com- 
ponents of the magnetic and electric vectors resolved 
in the same direction, that of the common axis of 
its two antennas. The mechanical design and lay- 
out of the satellite are shown in figure 4. 




Figure 4. The mechanical design of the satellite {the outer 
cover is not shown) . 



In order to explore the structure of the wave field, 
the satellite would be set spinning, so that the direc- 
tion of the ai it oinui axis, and hence the amplitude 
and phase of the received signals, would change 
continually. Information is required on the vari- 
ations of the amplitudes of the signals and of the 
sign of their correlation. These variations would 
be easiest to interpret if the spin axis of the satellite 
was at right angles to the axis of the antennas; then 
the motion of the antenna axis would be a simple 
rotation in a plane, and all variations would be 
periodic. The theory for interpreting these vari- 
ations, as developed in sections 4 and 5, assumes 
that this relation between the spin axis and the 
antenna axis is established exactly. 

Details of various aspects of the design of the 
satellite now follow. 

3.2. VLF Antennas 

The design of the loop antenna is quite straight- 
forward. At vlf, the most efficient arrangement is 
an air-cored loop, connected to the receiver through 
a matching transformer. The basis for optimum 
design, subject to limitations of weight and space, 
is provided by the work of Helliwell [4]. 

The design of the electric dipole is more difficult, 
since this antenna is in direct contact with the iono- 
spheric plasma, and its properties are affected 
thereby. The aims of design are to ensure that the 
effective height of the antenna is the same in the 
medium as in free space, and that its impedance 
does not vary excessively as the satellite spins. 
Consideration of the probable effects of the medium 
suggest that the best arrangement is one where each 
arm of the dipole consists of a conducting sphere 
supported by a long insulating rod; connection is 
made to the sphere by a thin wire running down the 
center of the rod. The merits of this arrangement 
are that the sensitive parts of the antenna are well 
localized, and are removed as far as possible from 
the disturbing influence of the body of the satellite. 
Furthermore, the effective height of this antenna 
should not be altered by the formation of a "positive- 
ion sheath" around the satellite (Seddon [5]), though 
its impedance doubtless would be altered. To mini- 
mize the effects of changes of antenna impedance, 
the input impedance of the receiver should be as high 
as possible. 

3.3. Electronic Equipment 

A block diagram of the electronic equipment is 
given in figure 5. It comprises vlf receivers, detec- 
tors, a telemetry system, and a generator for cali- 
bration signals. 

Two separate vlf receivers are needed, one for 
each antenna. They should be designed for low 
noise, stability of phase shift, and wide dynamic 
range. Their pass-bands should be fairly narrow, 
to reduce atmospheric noise, but not so narrow that 
the phases of their outputs fail to reproduce the 
changes in phase of the received signals that accom- 
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Figure ">. The electronic equipment of the satellite. 



pany the spinning of tln^ satellite. They should 
incorporate negative feedback to improve both the 

stability and the dynamic range. 

The amplitude of the signal that emerges from 
each receiver is measured by a detector, the output 
from which is telemetered to ground. Also, the two 
signals are compared in a phase detector, and the 
polarity of the output relayed to ground on a third 
telemetry channel. The ground station records 
these three data, which are proportional to the 
amplitudes of the received signals and the sign of 
their correlation. The factors of proportionality 
would have to be determined by calibration of the 
whole system prior to flight. Also, to confirm that 
the system is working properly during flight, it 
should be calibrated automatically from time to 
time by injecting standard signals at the inputs to 
the receivers. 

3.4. Orbital and Spin Motion 

The orbit of the satellite should be eccentric, so 
that measurements are obtained over a range of 
height; an apogee height of 1 to 2000 km and a 
perigee height of about 300 km would be suitable. 
At this height, the perigee would lie below the maxi- 
mum of the F2 layer from time to time, and then 
the measurement of electron density could be 
checked by comparison with ionospheric soundings. 
There are no special requirements for the inclination, 
except that it must be large enough to allow the 
orbit to pass over the transmitter, and that the value 
63.4°, at which the line of apsides does not precess, 
should be avoided. 

The necessary rate of spin is determined by tin' 
requirement that the satellite, in moving along its 
orbit, should not experience much change in local 
electron density during one rotation. If the dis- 
tribution of electrons in the outer ionosphere is 
assumed to contain no irregularities smaller in scale 



than 1 km, and the speed of the satellite is taken to 
be roughly 5 kmps, then the rate of spin should not 
be less than 5 rps. If the smallest irregularities 
actually are Larger than assumed here, then propor- 
tionately Lower rates would be acceptable. 

The absolute direction of the spin axis in space 
should be chosen so that the satellite is oriented 
favorably for making measurements at the times 
when it passes over the transmitter. At such times, 
the spin axis should be more or less horizontal, and 
should lie roughly in the plane of the orbit. It has 
to be horizontal so that the spinning motion will 
bring the antenna axis through the vertical, where 
the electric dipole is coupled strongly to the field 
component h\: for it is on this coupling that the 
measurement of the angle d depends (see sees. 4 and 
5). The spin axis also has to lie close to the plane 
of the orbit, to ensure that the elements of the elec- 
tric dipole do not pass behind the body of the satel- 
lite 4 , where they might encounter a wake of disturbed 
ionization. However, it is doubtful whether such a 
wake would be formed in the outer ionosphere, where 
the mean free path of the particles is much larger 
than the dimensions of the satellite, so this require- 
ment may not be crucial. 

3.5. Mechanical Design 

Bracewell and Garriott [6| have considered the 
influence of mechanical design on the free rotation 
of a satellite. They point out that the character- 
istics of the i'v^v rotation are determined by the 
relative values of the principal moments of inertia, 
and that if one moment is much larger than the other 
two, then the satellite will tend to spin about the 
axis associated with that moment. Therefore, the 
mass of the satellite should be arranged so that its 
moment of inertia is greatest about the desired axis 
of spin. If a satellite of this type was set spinning 
initially in a direction slightly incorrect, then its 
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subsequent motion would include a wobble at twice 
the spin frequency. This wobble could be damped 
out by such means as suitably disposed tubes of 
viscous fluid, whereupon the satellite would be left 
spinning smoothly about the correct axis. 

In the present case, the main aim is to ensure 
that the spin axis is transverse to the axis of the 
vlf antenna system. If this were the only aim, and 
the particular location of the spin axis in the trans- 
verse plane could be arbitrary, then it would follow 
that the satellite should be long and thin, extended 
along the axis of the vlf antenna system. However, 
it is better to fix the direction of the spin axis in the 
frame of the satellite completely, so as to provide 
a stable direction for the antenna of the telemetry 
system; to this end the mass should be distributed 
as in a disk, to make the moment of inertia greatest 
about the chosen direction. 

Once the spinning motion has been established 
about the desired axis in the frame of the satellite, 
the absolute direction of this axis in space should 
remain fixed by gyroscopic action. However, if any 
force were to exert a systematic couple on the satellite 
as it moved along its orbit, then the spin axis would 
precess. Possibly the forces due to atmospheric 
drag, or to electric currents induced in the body of 
the satellite as it spins in the earth's magnetic field, 
could exert such couples. The couple due to atmos- 
pheric drag could be eliminated by designing the 
satellite with reflection symmetry about a plane 
perpendicular to the spin axis, while induced currents 
could be prevented by making the body of the 
satellite in several parts that were insulated from 
one another. 

A mechanical design that embodies these principles 
is illustrated in figure 4; in this drawing, the protec- 
tive outer shell of the satellite is supposed to have 
been removed. 



4. Received Signals and Their Variations 
4.1. General Behaviour 

This section deals with the signals received on the 
two vlf antennas, and how they vary as the satellite 
spins. Specifically, expressions are derived for the 
amplitudes of the two signals and for the sign of 
their correlation, since this is the information that is 
relayed from the satellite to the ground. 

Now the loop antenna and the electric dipole 
antenna are set on a common axis (sec. 3), so the 
signals that they respond to are, respectively, the 
components of H and of E along this axis; call these 
components H a and E a . Take the coordinate system 
of figure 2, and let the instantaneous direction of the 
antenna axis be specified by its direction cosines 
&!, a 2 , and a 3 relative to the coordinate axes. Then 
the received signals are 



magnetic signal: H a —aiHi-\-a 2 H 2 
electric signal: E a =a l E 1 -{-a 2 E 2 +a 3 E 3 . 



(24) 
(25) 



It is convenient to express H a in terms of H 2 alone, 
and E a in terms of E u while recalling that the field 
components H 2 and E x are in phase, and that their 
ratio is the wave admittance (sec. 2.4.). In elim- 
inating the other components, the structure of the 
field will be assumed to be that given by the approxi- 
mations of section 2.5. After making the appro- 
priate substitutions, the following expressions are 
obtained for the received signals: 

magnetic signal: H a = {a 2 ±ja x }H 2 (26) 

electric signal: E a ={a x — a z t&ndTja^Ei. (27) 

The sign depends on the direction of propagation 
with respect to the geomagnetic field; the upper sign 
applies when there is a component of propagation 
along the field. 

As the satellite moves along its orbit, and spins 
about its axis, the received signals vary in amplitude 
and phase. The variations due to the orbital motion 
are slow and systematic, while those due to the 
spinning motion are relatively rapid and are, for 
the most part, periodic. In considering the varia- 
tions over a single rotation period of the satellite, 
only the effects of spin need to be considered. 

One general feature of these variations can be seen 
immediately; it is that they all must be periodic at 
twice the spin frequency, since a reversal of the 
direction of the antenna does not alter the ampli- 
tudes of the received signals, nor their relative phase 
(except possibly by 360°). 

To find the exact form of the variations, expres- 
sions must be obtained for the direction cosines of 
the antenna axis as functions of time. 



4.2. Direction of the Antenna Axis 

The geometry of the situation is represented in 
figure 6. Figure 6a is a perspective drawing that 
gives the notation for the angles that specify the 
directions of the various vectors. In figure 6b the 
same information is given in the form of a stereo- 
graphic projection onto the plane of the wavefront. 
Details of this projection can be found in any stand- 
ard work on crystallography, such as that of Tutton 
[7] or Bunn [8]. It maps the surface of a sphere into 
the interior of a circle on a plane, and has the prop- 
erty that arcs of great circles on the sphere are 
represented by arc of circles on the plane; those 
representing great circular arcs on the upper hemi- 
sphere are drawn as solid lines, while the projections 
of arcs on the lower hemisphere are drawn as broken 
lines. Points in the upper hemisphere are repre- 
sented by solid dots, and in the lower hemisphere by 
small open circles. In the present application, the 
orientations of vectors and planes in space are repre- 
sented by the projections of their intersections with 
an imaginary sphere that is concentric with the 
satellite. 
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Figure 6. The geometry of the spinning motion. 
(a) In perspective, (b) on a stereographic projection. 



The direction of the spin axis of the satellite is 
taken as that direction about which the rotation is 
right-handed. It is assumed to be perpendicular 
to the antenna axis, as specified in section 3. The 
angle that it makes with axis 3 of the coordinate 
system is called SF, and the angular velocity of the 
spinning motion is called 12. 

The time t is measured from an instant when the 
antenna axis lies in the plane of the wavefront 
(the 1-2 plane). The position of the antenna axis 
in this plane is then given by the angle 0, which is 
measured from axis 1 to that end of the antenna 
axis that is going upwards, toward axis 3, at this 
particular instant. It is reckoned positive if it is 
measured toward axis 2, and negative otherwise; 
thus, in figure 6, </> is positive. 



In these terms, the three direction cosines are 
given by the expressions 



ai = coS(f) costtt— sin</> cos\I> sinQt 
a 2 =sincl) cosf#+cos0 cos^ sinilt 
a 3 =siml>' sinS2£. 



(28) 
(29) 
(30) 



The angles ty, </>, and 6 define completely the rela- 
tive directions of the spin axis, the wave normal, 
and the geomagnetic field. They vary progressively 
as the satellite moves along its orbit, but these 
variations are slow compared to the spinning 
motion itself. 

4.3. Amplitude of the Magnetic Signal 

The real and imaginary parts of the magnetic 
signal, relative to II 2 as the standard of phase, may 
be found by substituting the expressions for the 
direction cosines into (26). Their values, normalized 
with respect to Z/ 2 , are 

R{H a /H 2 ] =a 2 =sin0 cosl#+cos0 cos^ sinQg (31) 

I{HJH 2 } = ±ai= ± (cos0 cosl#— sin<£ cos^ sinQt) 

(32) 

Thus the square of the normalized peak amplitude of 
tta is 



\HJH^=a\+aS 

= 1 — sin 2 ^ sin 2 fi£ 



(33) 

(34) 



Note that this expression does not involve either <t> 
or 0. It is made up of a constant part, and a part 
that varies sinusoidally at twice the spin frequency. 
The normalized peak amplitude, given by the 
square root of (34), takes its maximum value of 
unity when t lie antenna lies in the plane of the wave- 
front (ttt=0, 7r), and its minimum value of |cos^| 
at points halfway in between (£#=7r/2, St/2). 

4.4. Amplitude of the Electric Signal 

Similarly, the real and imaginary parts of the 
electric signal E a , normalized with respect to E u are 



R { EalE x }=a 1 —a 3 tan0 

= cos0 coslW— (sin0 eos^+tan0 sin^) 
sinl# 



(35) 



I{£ , a /£ , 1 } = : Fa2 ==: F(sin4> cosQt-\-cos<f> cos* siniM)- 

(36) 

Hence the square of the normalized peak amplitude 

is 



|# 8 /J2i| t =(a 1 -a,tan0)*+a* 2 



(37) 



= 1 — (1— tan 2 0) sin 2 * sin 2 fi< — 2 tan0 sin* 

sinQ<(cos</> cosilt— sin</> cos* s'mQt. (38) 
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This is more complicated than the corresponding 
expression for the magnetic signal, because the elec- 
tric field has a longitudinal component. However, 
it also represents the sum of a constant term, and a 
sinusoidal term of angular frequency 212. The nor- 
malized amplitude of E a also is equal to unity when 
Qt=0 or t, while its value at the points halfway in 
between is 

|S a /Ei| 2 =cos 2 ¥+tan0 sin0sin2^+tan 2 0sin 2 ^. (39) 

In general, however, these are not its maximum and 
minimum values. The times of stationary amplitude 
are separated by one quarter of a rotation, and are 
given by 

= 1 cot -1 {sec0 (sin<£ cos^— cot20 sin^)}. (40) 

The question of which time corresponds to the maxi- 
mum may be decided by inspecting the geometry of 
the situation. In most situations, maximum ampli- 
tude occurs near to the time when the antenna axis 
passes through the plane perpendicular to the 
geomagnetic field, because this is the plane of the 
electric vector. 

4.5. Sign of the Correlation 

Here the term "correlation" is being used rather 
loosely, to mean the average of the instantaneous 
product of the signals H a and E a , taken over one 
complete cycle of the radiofrequency. If the two 
signals were in phase, its value would be ^ \H a \ \E a \. 
In general, it is 



H a E a =m{H a }R(E a }+I{H a }I{E a }] (41) 



:f \H a \ \E a \ COS (a H — OL E ) 



(42) 



where a H and a E are the phases of the magnetic and 
electric signals, respectively, measured in relation to 
the phase of H 2 or E x . From (42), it appears that 
the correlation is positive when (a H —a E )<C90°, nega- 
tive when {a H — a^)>90°, and zero when the two 
signals are exactly in phase quadrature. 

The correlation could be determined by applying 
the voltages that represent H a and E a to the inputs 
of an electronic multiplier, and then smoothing the 
output with a time constant that is long compared 
to the period of the radiofrequency, but short com- 
pared to the period of the satellite spin. The pro- 
posed satellite would probably use some simpler 
type of phase detector, rather than a true multiplier. 
The output from such a detector would not neces- 
sarily be proportional in magnitude to the correla- 
tion, but it would have the same sign. For the 
present purpose, only the sign is of interest, so it is 
satisfactory to discuss the correlation. 

The correlation of the received signals is given as 
a function of the direction cosines by the expression 



H a En 



which contains the correlation of H 2 and E x as a 
normalizing factor. From this expression it appears 
that the sign of the correlation depends jointly on 
two circumstances: first, whether 6 is greater or less 
than 90°, and second, which of four quadrants of the 
sphere contains the direction of the antenna axis. 
These quadrants are bounded by the 1-2 plane and 
the 1-3 plane (fig. 7a). When the antenna axis lies 
in either of these planes, passing from one quadrant 
to another, the correlation is zero. The rules for 
the signs are given in figure 7b. 




-*2 



U90° 



?>90° 



Figure 



■\H t \ \E\ 



= —a 2 a 3 tan0 



(43) 



How the sign of the correlation is governed by the 
direction of the antenna axis. 

(a) The four quadrants of the sphere, within each of which the sign is constant 
(b) the rules for the sign. 

These results can be explained as follows: The 
departures of the relative phase of H a and E a from 
the quadrature condition are due entirely to the 
longitudinal component E 3 of the electric field. If 
this component was absent, as indeed is the case 
when 0=0, then the total E vector would be wholly 
transverse and would execute the same motion as 
the H vector, with a phase lag or lead of 90°; the 
received signals H a and E a , would then be in phase 
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quadrature for all directions of the antenna axis. 
When E\ is finite (05^0), the signals are exactly in 
quadrature in just two special situations: the first is 
the situation where the antenna axis lies in the 1-2 
plane (1^=0, ir), so that the electric dipole is not 
coupled to E 3 : the second is the situation where the 
antenna axis lies in the 1-3 plane, so that the dipole 
is not coupled to E 2 . Then, since E z is in phase with 
E x (see sec. 2.5), the electric signal has the same 
phase as it would if E 3 were absent. The times of 
passage of the antenna axis through the 1-3 plane 
are given by 

f2/ = tan" ^-sec^ tan0}. (44) 



4.6. An Example of the Variations 

The example of figure 8 shows how the properties 
of the received signals are expected to vary during 
one half-revolution of the satellite, for a particular 
orientation of the spin axis. Figure 8a is a graph 
of the normalized values of the signal strengths 
\H a \ 2 and \E a \ 2 , and of the correlation H a E a , as 
observed at the satellite; note that these three quan- 
tities all vary sinusoidally. Figure 8b shows how 
th e corr esponding records of \H a \, \E a \ } and the sign 
of H a E a would appear at the ground station. 





Figure 8. The variations of the properties of the received 
signals, for the particular case where ^ = 40°, <t> = 30°, and 
6 = 160°. 
(a) The variations of |flal 2 , \E a \ 2 , and E a H a , (b) the corresponding records. 



5. Analysis of the Records 

5.1. Wave Admittance, Refractive Index, and 
Electron Density 

Records such as that of figure 8b contain almost 
all the information that is needed for calculating the 
electron density of the ionosphere in the neighbour- 
hood of the satellite. The calculation, though, is 
somewhat indirect. 

The quantity that can be obtained most directly 
from the records is the local wave admittance, A. 
This, it will be recalled, is the ratio of the field com- 
ponents H 2 and E x . However, since the polarization 
in the plane of the wavefront is circular, both for H 
and for E (sec. 2.3), the wave admittance is uiven 
also by the ratio of the electric signal E a , taken at 
any instant when the antenna axis passes through 
this plane. Such instants can be recognized from 
the fact that the magnetic signal is then at a maxi- 
mum. 

Given the wave admittance, the local refractive 
index /jl can be calculated immediately from (2:3). 
Thus 

n=A/A . (45) 
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The relationship between ju and the local electron 
density N is contained in (19) . After rearrangement 
to express N as a function of the other variables, this 
equation reads 



AT = (^).|h |^ C os^co}(m 2 -1). 



(46) 



A convenient approximation to this expression may 
be obtained rearranging (6) similarly; it is 



N= 



H \cos 0|ju 2 



(47) 



in which c is the speed of light. Clearly, to calculate 
the electron density from a given value of fx, it is 
necessary to know also the local magnetic field 
strength H , and the angle between the direction 
of propagation and the field. 

The magnetic field strength cannot be determined 
from the records. If the position of the satellite is 
known, the field strength at this point can be esti- 
mated roughly from the known distribution of field 
on the earth's surface, by extrapolating upwards. 
It would be better, however, to measure the field 
directly by means of a magnetometer borne in the 
satellite; this possibility is discussed in section 6. 

The angle 6 can be obtained from the details of the 
variations of the recorded data. The argument, 
which is rather circuitous, depends on the following 
results from the previous section: (a) The variations 
of the amplitude of the magnetic signal are governed 
by SI> only (34), (b) the times of change of the sign 
of the correlation are governed by \I> and <£, but are 
independent of 6 (44), (c) the variations of the 
amplitude of the electric signal depend on all three 
angles, including 6 (38). 

Evidently, by analysing these three variations in 
the given sequence, it is possible first to calculate the 
angle ty, then the angle <£, and finally to obtain the 
angle 6. These calulations are discussed in the rest 
of this section. The discussion is illustrated by 
analysing the record of figure 8b. 

5.2. The Angle ^ 

\I> is the angle between the spin axis and the wave 
normal (sec. 4.2), and its range of possible values is 
to 180°. The value of ^ may be obtained from the 
observed variations of the amplitude of the magnetic 
signal ; in fact, from (34), it is given by the ratio of the 
extreme values of the amplitude. That is, 



|cosSF| = 



\H a \ min 



\H a \ max 



(48) 



Hence ^f is determined, though with ambiguity. For 
a given value of |cos ^|, there are two possible values 
of ^; their sum is 180°. 

Thus, in the example of figure 8, the ratio of the 
minimum and maximum values of |fiT a | is 0.766; on 
this evidence, ^ is 40° or 140°. The possible loca- 
tions of the direction of the spin axis are shown on 
the stereographic projection of figure 9. 




Figure 9. The possible directions of the spin axis, for a given 
value of |cos^|. 

5.3. The Angle 4> 

<t> is measured from axis 1 to the upgoing end of the 
antenna axis, at an instant when the latter passes 
through the 1-2 plane (sec. 4.2); its possible values 
range from —180° to +180°. The value of <£ may 
be obtained from the times of change of the sign of 
the correlation. Recall, from section 4.5, that the 
correlation changes sign at two instants in each half- 
revolution of the satellite: the first, which can be 
identified readily since it is also the instant when 
|i7 a | is a maximum, occurs when the antenna axis 
passes through the 1-2 plane (tit=0): the second 
occurs when it passes through the 1-3 plane. If the 
observed time delay between the two instants is 
expressed, in terms of the period of rotation of the 
satellite, as a phase angle tit, then <j> can be calculated 
by using (44). After rearrangement, this equation 
gives 

tan $= — cos ^ tan tit. (49) 

Altogether, there are now four possibilities, since 
for each value of ty there are two values of <j> that 
satisfy (49), and \[/ itself has two possible values. 
The problem now is how to resolve this ambiguity 
and select the correct combination of \f/ and <t>. 

5.4. Resolution of the Ambiguity 

If their true values are denoted by the symbols 
M>* and <j> t , the four combinations of ^ and <j> that 
would be possible in the light of the presentjevidence 
are the following : 

(a) \p t <$>t 

(b) + t 0-180° 

(c) 180°-*, 18O°-0 f 

(d) 180°-*, -<t> t . 
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Figure 10 illustrates these four possibilities in the 
analysis of the record of figure 8b; figure 10a shows 
the possible directions of the spin axis, and figure 
101) the corresponding loci of the antenna axis. 

What allows this ambiguity to be resolved is the 
fact that the angle 6 is not completely unknown at 
the outset of the experiment. Thus, in the example, 
6 is certainly greater than 90°, because the trans- 
mitter is located in the northern hemisphere where 





Figure 10. The fourfold ambiguity in the determination of the 
direction of the spin axis. 

(a) The possible directions of the spin axis, (b) the corresponding loci of the 
direction of the antenna axis. 



the lines of magnetic force point downwards, while 
the waves are traveling upwards. Granted this 
knowledge, the alternatives (b) and (d) can be 
eliminated by considering the direction in which the 
correlation changes sign at ilt = (); in the example, 
the sign changes from negative to positive at this 
instant . Now, since 8 is greater than 00°, the sign is 
governed by the rule on the right -hand side of 
figure 7h. Applying this rule, it appears that ;it 
the instant 12£ = the antenna axis is entering the 
upper right-hand quadrant of the sphere, and hence 
its locus must be either (a) or (c). 

The angle can be estimated more closely by 
assuming that the wave normal is directed to the 
vertical by refraction (see sec. 2.1), and by making 
use of the known value of the magnetic dip; at a 
guess, such an estimate would lie witliin 10° of the 
true value. If so, then the choice between the two 
remaining alternatives can be made by considering 
the variations of the amplitude of the electric signal. 
The expected variations of 12^ | are calculated for the 
two possible loci, using the estimated value of 0, and 
they are compared with the observed variation. 
The two calculated variations are mutually opposite. 
so that usually one should fit the record much better 
than the other, in spite of the inaccuracy of the 
estimate of 6. Figure 11 shows the comparison for 
the given example, taking 150° as the estimate of Y. 
This estimate is in error by 10°, but nevertheless 
it is clear that (a) is the correct locus. Thus, the 
ambiguity is resolved. 

In practice, the direction of the spin axis in space 4 
should be fairly well known from the circumstances 
of the launching, so that the ambiguity could be 
resolved immediately, and the chain of reasoning 
described above would he superfluous. 

5.5. The Angle 6 

6 is the inclination of the wave normal to the 
direction of the earth's field, and its possible values 
lies between 0° and 180°. Now that ^ and are 
known, 6 may be found by adjusting its value until 
the observed and calculated variations of the ampli- 
tude of the electric signal agree most closely. 

When 6 is known, the local electron density N 
can be calculated from the equations given in 
section 5.1. 

Since 6 is estimated from the details of the modula- 
tion of the received signals by the spinning motion 
of the satellite, obviously the depth of this modula- 
tion is one of the main factors that control the 
accuracy of the estimate. Now the depth of the 
modulation is governed chiefly by the inclination of 
the spin axis to the wave normal, that is to say, 
by the value of M>. The most favorable condition 
is that ^^90°, which gives 100 percent modulation 
of the magnetic signal; in section 3.4, the absolute 
direction of the spin axis in space was chosen so as 
to fulfill this condition when the satellite passes over 
the transmitter. On the other hand, the least favor- 
able condition is that ^^0°, where neither the 
magnetic nor the electric signal is modulated. If 
this latter condition were to arise, then 6 could not 
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Figure 11. Comparison of the "observed" and calculated variations of the strength 
of the electric signal: in the calculations, it was assumed that = 150° . 



be determined from the records. A theoretical esti- 
mate of 0, obtained as described in section 5.4 above, 
would then have to be used in calculating the electron 
density N. Note therefore, from (37), that the 
inferred value of A T is roughly proportional to |cos0|, 
so that it should be insensitive to errors in d around 
the assumed value of 160°. 

6. Integrated VLF Experiment 

In the experiment, as proposed so far, the satellite 
measures the least number of quantities that suffice 
to determine the electron density. Much more 
could be learned about the ionosphere if, at the same 
time, measurements were made of certain other 
quantities; namely, (a) The aspect of the satellite, 
(b) the local magnetic field, (c) a component of the 
vlf wave field (H or E) parallel to the spin axis. 
These measurements will be discussed in turn. 

(a) Aspect. The instantaneous aspect of the satel- 
lite could be measured by photoelectric devices 
which sense the direction of the sun, the moon, or 
the horizon of the earth. Such measurements would 
determine the absolute direction of the spin axis 
in space. 

This information, when combined with the record 
of the amplitude of the magnetic signal H a} would 
determine the absolute direction of the wave-normal. 

(b) Magnetic field . Three components of the local 
magnetic field could be measured by a pair of 
magnetometers. One, a fluxgate, would measure the 
component of field along the spin axis. The other, 
an induction magnetometer, would be set up in the 
plane transverse to the spin axis; it would make use 
of the spin of the satellite to induce an alternating 
voltage that would represent the two components 
of the magnetic field in the transverse plane. A 
magnetometer of this type has been used successfully 
in the lunar probe "Pioneer I" (Sonnett [9]). 



From the measurements of the three components, 
the total local field strength (H ) can be calculated. 
These same measurements, when combined with the 
measurement of aspect, also determine the absolute 
direction of the field. 

Having determined the absolute directions both 
of the field and of the wave-normal, the angle (6) 
between them can be obtained directly. This value 
is likely to be more accurate than that calculated by 
the method of section 5, where much use was made 
of the sign of the correlation between H a and E a ; this 
measurement is apt to be sensitive to errors caused 
by small phase-shifts in the receivers, and by noise 
also. 

It will be recalled that the values of the H and 6 
are needed for calculating the local electron density 
from the measured wave admittance (sec. 5.1). 

(c) Component of vlf field along the spin axis. 
Finally, following a proposal by Helliwell [10], the 
satellite could be equipped with a third vlf receiver, 
connected to an antenna that was alined with the 
spin axis. The signal received on such an antenna 
would not be affected by the spinning motion. Its 
phase (relative to a stable local standard carried in 
the satellite), and also its amplitude, would be tele- 
metered back to the ground. 

The rate of change of phase, when corrected for 
the change in path length for the telemetry trans- 
mission, yields the Doppler frequency shift of the 
vlf signals as observed at the satellite. This fre- 
quency shift is given by 



A/=(VCosS)/\ 



(50) 



where V is the speed of the satellite along its orbit, 
<5 is the angle between the orbit and the wave-normal, 
and X is the wavelength in the medium. Now the 
orbit is presumed known, and the local direction of 
the wave normal is measured, so the values of V 
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and <5 are both data. Hence, if A/ is appreciable, X 
can be calculated from its observed value. Now 



X = X //x 



(51) 



where \ is the wavelength in free spaee. Given X 
and \o, fJL can be calculated. Thus the Doppler 
shift would provide another measure of the local 
refractive index, and hence of the electron density; 
this measurement would check the value calculated 
from the wave admittance (sec. 5.1). 

If the vlf transmissions were modulated in some 
way, the modulation would be received at the 
satellite with an appreciable delay, because the group 
velocity of the whistler mode is rather low. The 
delay could be measured most easily on this third vlf 
channel, where the received signals are not modu- 
lated in addition by the spin of the satellite. From 
the observed delay, it is possible to determine the 
integral of the plasma frequency with respect (o 
height up to the level of the satellite (Ilelliwell [10]). 

The proposal is made that there should be an 
integrated vlf satellite experiment, incorporating 
all the features described above. The following 
information about the ionosphere would be ob- 
tained at each point on the orbit: (a) Local strength 
of earth's magnetic held, (b) local direction of field, 
(c) local strength of vlf waves, (d) local direction of 
propagation, (e) local electron density, (f) height 
integral of plasma frequency. 

The changes in these quantities, as the satellite 
moves along its orbit, would be expected to show 
both systematic and random components. The 
systematic components would be related to the bulk 
properties of the ionosphere, while the random 
variations would give information on ionospheric 
irregularities. 

7. Discussion 

The main proposal made in this paper is for a 
satellite that would determine local electron density 
in the ionosphere, by exploring the structure of a 
vlf wave field set up by a transmitter on the ground. 
The following sources of difficulty are anticipated: 

(a) Noise. One limit to the accuracy of the ex- 
periment will be set by noise. Judging by previous 
experience at vlf, there should be no difficulty in 
making the receivers sufficiently quiet that the limit 
is set by the external atmospheric noise picked up 
on the antennas, which is unavoidable. However, 
the actual level of atmospheric noise within the 
ionosphere is hard to predict. At the ground, vlf 
noise is observed to come both from the lower at- 
mosphere, where it is produced by lightning strokes, 
and also from the ionosphere, where its mode of 
origin is uncertain (Watts [11|). Within the iono- 
sphere, presumably, the former source of noise is less 
effective and the latter more so. It is hoped that 
the vlf satellite experiment already planned by 
Helliwell [10] will provide data on the noise level in 
the ionosphere. 



(b) Multiple waves. The theory has assumed 
that the vlf field in the ionosphere is that of a single 
plane wave (sec. 2.1). Thus the experiment would 
not be expected to work near the conjugate point, 
where the downcoming direct waves would be 
mixed with upgoing waves reflected from the base 
of the ionosphere. Difficulty might be experienced 
even near the transmitter, if the downcoming wave 
returned from the opposite hemisphere was at all 
comparable in strength to the upgoing primary 
wave. This difficulty could be overcome either by 
working at a higher frequency, which would be more 
heavily absorbed, or by modulating the transmis- 
sions in such a way that the direct signal and the 
echo never overlap one another; for 15.5 kc signals 
from the station NSS, the time of travel over the 
whistler route and back would be about 1 sec (Helli- 
well and Gehrels [12]), so that a suitable form of 
modulation would he a % sec pulse repeated every 
2 sec. A potentially more serious problem is pre- 
sented by ionospheric irregularities, which would 
distort the upgoing wavefronts. If the curvature 
of the wavefronts was appreciable in terms of a 
wavelength, then the plane-wave theory would not 
apply. The significance of the measurement of 
wave admittance in these circumstances has not 
yet been investigated. 

(c) Effect of the medium on the electric dipole 
antenna. There is doubt as to how well the electric 
dipole antenna will perform when immersed in the 
anisotropic plasma of the ionosphere. The design 
of antenna that was described in section 3.2. is in- 
tended to minimize the adverse effects of the medium. 
If such effects persist nevertheless, their presence 
could be detected from the fact that they will almost 
ceit a inly cause the variation of \E a \ 2 with time to 
depart from the pure sinusoid predicted by the theory 
(sec. 4.4). 

(d) Disturbance of the medium by the satellite. 
There is a chance that the satellite, by virtue of its 
rapid orbital motion, may disturb the medium ser- 
iously and thus alter the distribution of magnetic 
and electric fields in its immediate neighbourhood. 
The purely hydrodynamic disturbance is likely to 
be small, and confined to the region behind the satel- 
lite (sec. 3.2). However, Kraus and Watson [13] 
have pointed out that the satellite, in its passage 
through the ionosphere, acquires an electrostatic 
charge, and that the field of this charge should 
produce a hydromagnetic disturbance ahead of the 
satellite. It is uncertain how serious this effect 
will be. In this connection, it is worth noting that 
instruments have been developed to measure the 
state of charge of a satellite (Krassovsky [14]), and 
doubtless one could be arranged to control some 
discharging device so as to maintain a state of 
electrical neutrality, and thus eliminate the effect. 

(e) Incorrect spinning motion. Finally, there is 
the possibility of failure to establish the desired 
spinning motion: the satellite might be launched 
either with no spin, or spinning too slowly, or again 
the spin might be fast enough but accompanied by 
a spurious wobble that would take 4 some time to 
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suppress (see sees. 3.4 and 3.5). If the spin was 
non-existent or very slow, then the experiment 
would be spoiled completely. A wobble would not 
affect the measurement of wave admittance, since 
the spin would still bring* the antenna axis through 
the plane of the wavefront (see sec. 5.1), but it would 
spoil the measurement of 6. In such a case the elec- 
tron density could still be calculated by using a 
theoretical estimate of 6 (sec. 5.4), and this procedure 
would probably not entail much error (sec. 5.5), un- 
less 6 was varying widely and at random along the 
orbit due to the effects of ionospheric irregularities; 
this last condition would be obvious from the records 
of signal amplitude. The presence of a component of 
wobble in the spin motion could be detected most 
easily if the satellite carried some device for measuring 
aspect. 

In view of these several difficulties that may beset 
the determination of local electron density by the 
wave admittance method, it is clear that this experi- 
ment should not be performed in isolation, but only 
as part of an integrated vlf satellite experiment such 
as that proposed in section 6, where some of the 
measurements would provide cross-checks on others. 



The author is indebted for advice from C. G. 
Little, R. M. Gallet, and R. S. Lawrence of this 
Laboratory, from J. W. Warwick of the High 
Altitude Observatory of the University of Colorado, 
and from R. A. Helliwell of Stanford University. 
The proposal for an integrated vlf satellite experi- 
ment arose out of a discussion with Dr. Helliwell. 
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